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and  thus  the  tallied  vehicle  was  chosen  to  ir.rct  the  design  criteria.  The 
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ratio  of  25  and  was  very  stable  in  heading,  while  meeting  all  other  design 
requirements. 
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Preface 

Despite  very  raold  developments  In  aerodvnamlc  theory  and  the 
oreat  Improvements  In  the  useuqe  of  this  theory  throuch  the  develop- 
ment of  modern  computer  facilities,  much  of  the  dasinn  work  on  unusual 
configurations  Is  a trial  and  error  process  that  utilizes  approximate 
mathematical  models  of  the  conf Iquratlon  to  eld  the  dmsloner  In  the 
analysis  of  those  conf louret Ions.  The  followinn  thesis  demonstrates  the 
application  of  this  deslon  process  to  a specific  deslqn  request  from 
the  Air  Force  Avionics  Laboratory.  The  trial  and  error  process  would 
have  been  more  error  than  trial  were  It  not  for  the  help  of  Mr.  Russel 
Osborno,  an  ennlneer  for  the  Air  Force  Flloht  Dynamics  Laboratory.  I 
would  also  like  to  thank  my  thesis  advisor  for  his  patience  and  for 
provldlno  some  necessary  direction.  Additional  thanks  no  to  the  Air 
Force  Institute  of  Techno  I oqy  Model  Shop  for  their  natlence  In  build  inn 
and  frequently  rnpairlna  the  test  modols  nocessarv  for  this  Invest I nation. 
It  qoes  without  saylnn  that  I owe  much  to  my  wife,  Scotty,  for  Insulatlna 
me  from  the  outside  world  while  wrttlno  this  report. 
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Notation 
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Wing  span 

Horizontal  stab! 1 1 zar  span 

Vortical  stabilizer  span 

Vehicle  equllbrlum  drag  coefficient 

Vehicle  drag  change  with  ancle  of  attack 

Vehicle  equilibrium  lift  coefficient 

Chango  In  lift  with  pitch  rate 

Vehicle  lift  change  with  annlo  of  at  lack 

Change  in  lift  with  angle  of  attack  rate 

Horizontal  stabilizer  lift  change  with  aerie  of  attack 

Wing  lift  change  with  angle  of  attack 

Change  In  rolling  moment  with  sides! Ip 
Chango  In  rolling  moment  with  roll  rote 
Chango  In  rolllnq  moment  with  yaw  rate 
Vehicle  coefficient  of  moment 
Zero  lift  coefficient  of  momont 
Chango  It  moment  with  pitch  rote 
Change  In  ronvint  with  angle  of  attack 
Chango  In  romont  with  angle  of  attack  rate 
Chunco  In  yawing  momont  with  sideslip 
Chang 2 In  yawing  moment  with  roll  rate 
Chango  In  yawing  momont  with  yaw  rate 
Chango  In  si  deforce  with  sideslip 

vl 


Chanqo  In  sldeforce  with  sideslip  due  to  vertical  stabilizer 


C Chanqe  ^ si  deforce  with  roll  rate 

VP 

C„  Chanqe  In  side  force  with  yaw  rate 

Vr 

C Equilibrium  wolaht  coefficient 

% 

c Wlnq  chord 

c Horizontal  stabilizer  chord 

H 

cy  Vertical  stab! I Iter  chord 

g Acceleration  of  gravity 

t*  Nondtnenslonal  rolllnq  moment  of  Inertia 

x 

A 

I*  Nondlmentlonal  pitch  I no  moment  of  Inertia 

Y 

I'j  Nondlmentlonal  yawing  moment  of  Inertia 

It  Static  margin 

n 

lfa  Body  horizontal  length 

I Length  from  vehicle  center  of  gravity  to  horizontal  stabilizer 

aerodynamic  center 

S Wing  Area 

SH  Horizontal  stabl I Izor  area 

SQ  Body  cross-sectlor.^l  area  at  point  of  maximum  thickness 

S Vertical  stab H Izor  area 

v 

T Oscillation  period 

V Tall  volume  cooV t Iclent 
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w Vehicle  weight 

Zy  Vertical  dlstanco  from  vehicle  canter  of  gravity  to  aerodynamic 

centor  of  vortical  stabilizer 

AR  Aspect  ratio 

T Winn  dihedral 

Y#  Flight  path  angle 

C Damping  coefficient 


vll 


GAE/MC/740-t 


Abstract 

An  unpowtred,  ungulded  vehicle  to  ba  dap  toyed  from  a 4,75"  fla-a 
tuba  was  designed,  Tho  da* I on  cr I tar  I a wars: 

1)  six  cubic  Inch,  tvo  pound  payload 

2)  minimum  glide  spaad  of  90  knots 

3)  qllde  ratio  In  oxcass  of  fiva 

4)  stabla  with  particular  amphasls  on 
haad I ng  stab  1 1 1 ty 

Two  basic  configuration?  wara  cons I dared  In  the  analysis:  the  first  was 
a tailless  vohlcie,  and  the  second  was  a conventionally  tallied  vehicle. 
It  was  found  that  the  tail  loss  vahlcla  could  not  neat  the  stability  re* 
qulrements  and  thus  tho  tallied  vehicle  was  chosen  to  meat  the  design 
criteria.  The  analysis  procedure  consisted  of  a computer  model  to  analyze 
vahlcla  stab II I tv,  wind  tjnnol  tests  to  calculate  aerodynamic  deriva- 
tives, and  fllgh*  tests  to  demonstrate  performance.  The  finalized  vohlcle 
had  a folding  wing  and  a flexible  tail-boom  to  meet  the  packaging  re- 
quirement, The  resultlnq  vehicle  had  a glide  ratio  slightly  In  excess  of 
25  and  was  very  stable  In  heading,  while  meotlnq  all  other  design 
requirements. 
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I • Ini roduct I on 

The  Problow 

In  an  ora  of  nvar  chant}  I no  nor ; a I tactics  and  counter-tactics,  a 
constant  soared  tor  now  woanons  Is  walntalnod,  The  doslon  of  ono  such 
weapon  systen  Is  the  subject  of  this  thesis,  The  Air  Fere o Avionics  Lah- 
oratory  reouosted  tho  flesion  of  an  unpnwereo,  unouidod,  air  erne  vohlcle 
capable  of  carry I no  n sly  cubic  trch  payload  we  I oh  I no  two  pounds  at 
a mlnl'-u"  olldo  speed  of  50  Knot;,.  To  achieve  wax  I run  dispersal  capabil- 
ity, each  vehicle  oust  mlntnln  Its  deployment  no-wind  headlno,  and  rust 
b«  of  nlnlrur*  vo'ure  50  that  r<x>  i-jr.  n.r*pers  ray  pi»  carried  bv  tpe  car- 
rier vnhlclo.  To  ellr.lnato  the  need  for  dovolop,-«>nt  and  qualification 
of  a now  s*oro,  tho  p’-p^osal  reoilred  that  thn  vehicle  ho  able  to  he 
carried  In  a 4,75  Inch  flaro  tube,  thus  presentlno  a slonlf leant  pachao- 
Ino  problem,  |n  sircirv  the  prop  lee  presented  was  to  deslnn  a vehicle 
sub  |ect  to  the  follow  Inn  constraints: 

I)  uncovered  and  unnuided 

?)  six  cubic  inch,  two  pound  payload  capacity 

3)  minimum  nlido  speed  o*  ‘>0  knots 

4)  olldo  ratio  in  excess  of  five 

5)  stable  with  particular  emphasis  on  hoadinn  stability 

A)  minimum  sl.»e  and  oac^.vrnahlt?  In  a 4,75  inch  tube. 
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The  Analysis  Procedure 

As  in  any  design  study,  the  method  used  to  create  the  ftnal  result 
was  one  of  iteration  among  .ir  best  solutions  to  each  of  the  design  crite- 
ria. This  analytic  process  used  the  most  simple  tnethods  suitable  to  the 
complexity  of  the  problem  and  testing  methods  utilizing  the  best  avail- 
able equipment.  The  procedure  used  In  this  design  effort  was  divided  Into 
four  parts:  an  extensive  literature  search  to  find  the  most  suitable  con- 
figuration for  the  required  mission;  the  creation  of  a useable  mathemati- 
cal model  to  describe  the  chosen  configuration;  wind  tunnel  tests  to  de- 
termine the  required  aerodynamic  coefficients;  and  flight  tests  of  the 
various  test  vehicles. 

Iteration  among  each  of  the  latter  three  areas  was  accomplished  to 
enable  the  chosen  configurations  to  bo  optimized  for  the  mission.  This 
particular  type  of  design  process  allowed  the  configuration  to  be  modi- 
fied during  the  process  with  little  interruption  of  the  test  sequence. 

Limitations  and  Assumptions 

The  design  process  did  not  include  certain  areas  that  would  be  of 
Interest  In  further  studies  of  the  design.  These  areas  Include  deployment 
aerodynamics  and  stability,  and  a structural  analysis  designed  to  deter- 
mine the  construction  materials  which  prov'de  the  best  compromise  between 
weight,  strength,  and  flexibility.  These  two  areas  were  beyond  the  scope 
cf  this  study  due  to  limitations  in  the  available  facilities  and  time. 

The  major  assumptions  in  the  analytical  modeling  concerned  linearity 
of  the  aerodyr  imics  and  of  the  equations  of  motion  used  in  the  stability 
analysis.  8oth  of  the  assumptions  are  normally  made  in  the  preliminary 
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deslqn  phase,  and  are  nertfcularlv  valid  for  a nlider  deslqned  to  main- 
tain a symmetrical  equilibrium  fllqht  condition.  The  assumption  of  linear 
aerodynamics  made  the  stability  analysis  possible  since  the  stability 
derivatives  were  then  possible  to  evaluate  as  functions  of  the  planform 
and  equilibrium  fllqht  condition.  Further  assumptions  concern  I nq  each 
particular  conflagration  will  be  discussed  as  they  arisa. 

Order  of  Discussion 

The  report  first  discusses  the  choice  of  conf tquratlons  to  be  con- 
sidered, and  then  describes  the  entire  analytic  process  for  each  of  the 
chosen  conf lauratlons.  Althouqh  3! I of  the  phases  for  each  conf Iquratlon 
were  carried  out  simultaneously,  the  description  of  the  entire  analysis 
procedure  for  one  conf iquratlon  at  a time  leads  to  a more  understandable 
f -osentatlon.  Specific  conclusions  concern  I na  each  conf Iquratlon  are 
presented  at  the  end  of  the  discussion  of  that  particular  conf Iquratlon, 
while  qeneral  conclusions  and  recommendations  are  presented  at  the  end 
of  the  report. 
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II.  THF  CONTinURATION 


Configurations 

The  Initial  choice  of  configurations  to  bo  considered  was  Important 
because  It  provided  dlrecSlon  to  the  entire  design  process.  It  was  pos- 
sible to  rapidly  eliminate  some  configurations  from  further  study  because 
of  their  Inherent  limitations.  In  particular,  both  balloons  and  para- 
chutes could  not  meet  the  regulrement  for  a forward  velocity  of  50  knots 
and  were  therefore  eliminated  from  consideration.  All  parafoil  devices 
had  to  be  eliminated  due  to  glide  ratios  possible,  although  In  every 
other  respect  they  presented  an  excel  lent  solution.  It  was  possible  that 
a hybrid  vehicle,  oart  balloon  and  part  rigid  wing,  could  have  met  the 
design  reguirements,  but  storage  difficulties  for  such  a vehicle  elimin- 
ated It  from  consideration.  As  a result  of  the  above  eliminations,  only 
rigid  wing  vehicles  were  studied  further. 

Three  possible  configurations  were  evaluated  for  possible  study  and 
two  with  their  modifications  were  chosen  to  be  studied.  The  three  config- 
urations studied  were:  an  all-wing,  a wing-body,  and  a conventional 
tallied  vehicle. 

The  all-wing  vehicle  was  first  considered  to  represent  the  best  pos- 
sible solution  to  the  overall  problem.  After  studying  possible  plan- 
forms;  however.  It  was  found  that  an  excessively  large  vehicle  would  be 
necessary  to  provide  sufficient  useable  payload  volume.  Such  a large 
vehicle  presented  a formidable  packaging  problem,  and  the  folding  mech- 
anism was  excessively  complex.  For  these  reasons,  the  all-wing  design 
was  not  considered  further. 
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The  wing-body  was  considered  to  be  the  next  best  solution.  This 
vehicle  was  more  easily  packaged  than  the  all-wlno  vehicle,  had  the  poten- 
tial for  higher  performance  than  the  conventional  configuration  due  to 
lower  expected  drao,  and  provided  a less  complex  fold  I no  nroblem  than  the 
conventional  vehicle  and  was  thus  expected  to  be  more  reliable.  The 
critical  design  point  for  this  conf Iguratlon  was  the  stability  require- 
ment, This  point  was  expected  to  present  problems  since  a study  of  avail- 
able literature.  In  particular  reference  I,  Indicated  that  tailless 
vehicles  were  often  only  weakly  stable,  If  stable  at  all,  in  some  flight 
modes.  If  this  design  reoulrement  could  be  met,  then  this  configuration 
appeared  to  present  the  best  solution  to  the  design  problem. 

The  least  desirable  configuration  from  preliminary  estimates  was  the 
conventional  vehicle.  The  design  was  the  most  complex  mechanical ly  and 
thus  probably  the  most  unreliable  of  the  two  possible  configurations. 

It  did  have  certain  advantages;  however.  It  presented  the  most  stable 
configuration,  and  had  the  advantage  of  being  better  known  so  that  the 
design  process  could  be  better  analyzed. 

Configuration  Geometry 

After  the  two  basic  configurations  were  decided  upon,  the  basic 
planform  parameters  were  obtained  from  an  analysis  of  the  packaging  prob- 
lem. The  more  complex  packaging  problem  was  the  conventional  vehicle  and 
so  It  was  chosen  to  determine  the  gross  dimensions.  To  facilitate  the  con- 
struction and  testing  of  the  wind  tunnel  and  flight  tost  models,  and  to 
aid  In  the  comparison  of  the  two  configurations,  as  much  commonality  as 
was  possible  was  maintained  between  the  vehicles.  Afier  many  cut-and-try 
"paper  doll"  models,  the  final  gross  dimensions  shown  In  figures  I and  2 


were  chosen 
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Although  the  chosen  airfoil  section,  an  NACA  441?,  did  not  represent 
tho  opt*r.'.»m  solution  to  the  wlno-body  oroblem,  It  did  present  a ready  so- 
lution to  the  folding  problem  as  a result  of  its  nearly  flat  bottom  sur- 
face. A wine  foldlno  mechanism  was  designed  for  this  wing  by  Mr.  Lance 
Look  of  LOM  Company  under  contract  from  the  Air  Force  Flight  Dynamics  Lab- 
oratory. Through  use  of  a si mole  tension  spring  system,  ihe  wing  was 
folded  at  two  points  on  each  semi-span  as  shown  In  floures  1 and  2. 

A low  drag  wlno  tip  shape  was  developed  using  Information  from  ref- 
erence 2.  The  horizontal  cross  section  of  the  body  wts  chosen  to  be  sym- 
metrical with  the  point  of  maximum  thickness  located  at  the  body  midpoint. 
This  profile  was  used  rather  than  a more  streamlined  shapo  In  an  attempt 
to  ensure  that  the  body  center  of  pressure  was  behind  the  center  of  qrav- 
Ity  for  the  entire  static  maroln  ranoo  so  that  directional  stability 
would  not  be  a problem. 

The  one  remalnlno  parameter  to  be  established  was  the  wing  sweep. 

In  the  convents  i a t vehicle,  a wing  sweep  angle  of  0 degrees  was  chosen 
to  provide  maximum  lift  at  lower  angles  of  attack  and  thus  to  reduce  the 
drag.  In  the  wing-body,  wing  sweeo  was  a necessity  to  ensure  sufficient 
static  margin  range  to  provide  longitudinal  stability.  Wing  sweep  moved 
the  aerodynamic  center  aft  and  allowed  greater  latitude  In  payload  loca- 
tion, After  making  the  assumption  that  the  wing  aerodynamic  center  was  at 
the  ,25  mean  aerodynamic  chord  and  that  this  chord  was  located  at  the 
area  center  of  the  semi-span,  a qood  assumption  In  view  of  the  method  of 
reference  3 In  locating  this  chord,  it  was  found  that  a sweep  angle  of 
15  deqrees  would  provide  a static  margin  ran  ie  In  excess  of  .30  while 
presenting  the  least  possible  compromise  with  performance. 
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Ml.  THE  WiNQ-BOPY  DESIP*? 


Alrfoi  I Desjjnn 

All  flight  vehicles  must  satisfy  tha  following  thraa  requirements  to 

ba  able  to  achteva  stable,  trimmed  flight:  a positive  C^,  a negative 

C^,  and  a Cm"0  at  soma  angle  of  attack  greater  than  zero  (Ref  4).  In 

order  to  provide  this  capability,  a modification  of  the  basic  NACA  4412 

airfoil  section  had  to  ba  developed  since  a posltlv®!.*  cambered  airfoil 

could  noT  provide  the  required  positive  Cm  and  no  other  airfoil  section 

o 

could  be  found  to  provide  the  required  aerodynamic  coefficients  and  meet 
the  requirement  that  the  bottom  surface  be  flat  for  packaging  purposes. 
Several  solutions  to  this  problem  were  available  among  which  were  util- 
izing geometric  twist,  upward  deflected  flaps,  or  a reflexed  camber  line. 
The  simplest  solutions  wore  either  flaps  or  a reflexed  camber  line  since 
twist  would  require  construction  of  a new  wing,  while  the  other  solutions 
only  meant  sectioning  the  wing.  The  best  solution  from  the  packaging 
standpoint  was  the  reflexed  camber  line  so  this  solution  was  used. 

One  major  assumption  was  made  in  the  development  of  this  wing;  that 
creating  an  airfoil  section  with  sufficient  reflex  would  provide  suffi- 
cient accuracy,  ionorinq  all  three  dimensional  effects,  so  that  the  actual 
wing  camber  line  could  be  refined  In  the  wind  tunnel  without  reconstruc- 
tion. Resulting  wind  tunnel  tests  proved  the  validity  of  this  assumption. 

The  theory  of  thin  wing  sections  used  in  the  airfoil  design  (Ref  5) 
was  Iterated  over  several  profiles  until  the  deslrod  properties  were  ob- 
tained. The  nocossary  winn  fy,  obtained  from  a summation  of  moments  about 
the  aerodynamic  center  assuming  a static  margin  of  .20  and  neglecting  drag 

effects  on  the  moment  balance,  was  .1450.  This  yielded  a vehicle  (^*0 
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at  an  arbitrary  10*  angle  of  attack.  Making  tha  assumption  that  tha 
airfoil  could  be  represented  as  shown  below  In  figure  3* 


r 

c 

} 

o 

« 

o 

Fig.  3 ^eflexed  Camber  Airfoil 


and  that  line  A was  parabolic  In  shape,  lines  B and  C were  straight,  the 
equations  for  each  line  were  solved  by  matching  boundary  conditions  at 
points  (pt ,m( ) and  (p3,m3).  The  slopes  of  the  lines  at  these  same  points 
were  obtained  In  gereral  terms  as: 


dy  (2x«p1)(fn2P1-m1P2)  + w^Cp^-p^) 
<**A  Pja<P2-P!> 


dy  m -m^ 
dxB  P*"pl 


dX£  p2-l 


(1) 

(2) 
(3) 


After  utl 


where: 


I zing  Glauert*s  equation  for 


-H' 


dn  cos  nO  dO 


IS 


x » 1/2  (cos  0 ♦!) 


the  Fourier  coefficients: 


(4) 

(5) 
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the  Fourier  coefficients  nocotsary  to  calculate  C,  and  C were: 


. ?.pi 

• *Lv' 


co»“l(2p.-l)  ♦ — <cos*l<2p  -I > - cos“*(2p,-l )) 

* p -p  1 * 

¥t  1 


<« ,1* »-«  p 7(1-0  > ♦ m (p  p -p  2) 

♦ ■»  I I L.  L .1,  I i.  .L  (*-cos“l  (2p  -I ) ) ♦ (6) 

Y(W 


♦ 1 4 *(->ln(cot 


>l(2pl-m)J 


2 ~m  mm 

A • — — ^ sln(cos*M2p  •>!))+  -*-*($ ln(cos“l  (2p  -I  ))-slnlcos“U2p  -I ) ) ) 

1 *Lv  2 vi  1 1 


(l-p  )(»  0 -m  p ) m (p  p -pj) 

♦ ■■— — - - ■ I.  (-sln(cos“1(2p1-l ))) 


(7) 


» P.-*V9  « 

♦ -*r I — <-  - cos"l(2p  -I)  4 (sln(cQs”l(2o  -I )))) 

vvv  2 — » — 1 — 


.ili 


sln(2cos~i(2p,~i ))  rn^-n,  fs*n(2co.«~l  I2pj-I ) )-sln(2cos~l  (2p,-l ) ) > 


LV 


2 » 1 

Ypi 


ll-p  )<m  p -m  p )+m  (p  p -p  •)  (-sln(2cos"2(2p  -l)>) 
1 2 112  ' 12  l 1 


(8) 


p ?<p  -p  ) 

y\  ■>  *1 


(m  p -m  p ) r-s!n(cos-l(2p  -I))  jln(3cos-M2p  - 
2 112  1 1 


1 p -Hi  p > r- 

21  12 

t 2<p  -p  ) 

I K2  Ki  L 


• l))T 
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I 

I 
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r 


whore: 


Cj  • 2i|o-  lil_^£jand 


(9) 


Now  us!*?  <he  of  C|  and  Cm  In  farm*  of  that*  coefficient*  tor  an 

arbitrary  10*  eng I a of  attack  and  for  no  reflex  eng  la,  vhe  valuas  of  tha 
coefficient*  became: 

C|  • 1.8405 
C - -.1297 

nl 

while  reference  5 showed  values  of: 

Cj  ■ 1.5500 
Cm  - -.1000 

which  demonstrated  that  the  proposed  model  of  the  airfoil  section  wus 
accurate  enough  for  the  purposes  of  this  deslgr.  After  numerous  Iterations 
of  the  equations  to  arrive  at  a shape  that  y I aided  the  desired  properties, 
the  final  result  was  an  airfoil  section  with  a reflex  angle  of  13*  that 
started  the  reflex  at  the  .8  chord  point  as  Illustrated  In  figure  4. 


r 

(.4,. 12)  (.8,0.) 

Flq,  4 Reflexed  Airfoil 


Although  the  above  section  characteristics  could  have  been  extended 

to  the  finite  wing.  It  was  unnecessary  since  wind  tunnel  rosults  would 

yield  more  accurate  results  and  would  be  readv  for  the  stability  analysis. 
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Stab  I Mty  Analysis 

The  stability  analysis  rested  on  +ho  fact  that  the  tonoltudlnol 
stability  analysis  eo.Jd  jc  separated  from  the  lateral  analysts.  The 
equations  of  notion  used  in  the  analysts  wore  developed  In  reference  4 
and  the  assumptions  that  allow  the  longitudinal  notion  to  be  separated 
from  the  latoral  motion  ere  discussed  at  lenoth  there.  For  the  purposes 
of  the  following  analysis,  all  of  the  assumptions  were  valid. 

Longitudinal  Analysis.  Since  static  longitudinal  stability  could  bo 
normally  achieved  by  ensuring  that  C*.  Is  negative,  the  analysis  of  the 
dynamics  of  the  vehicle  also  Included  static  stability.  Tho  solution  to 
the  resulting  eigenvalue  problem  yloldort  nil  of  tho  roqulrod  Information 
concerning  dynamic  stability.  The  characteristic  matrix  resulting  from 
the  equations  of  motion  as  developed  In  roference  4 and  offer  eliminating 
both  thrust  and  speod  dor I vat I vu  terms  was; 


C slny 

% * 

Cl  “ co 

e o 

0 

-C«  cosy, 
e 

u 

— % — 
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a e 

2u  - CL 

d 

-C  Slny 

e * 
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The  assumption  that  tho  speed  derivatives  war*  negligible  was  warranted 
by  tha  vary  low  spaed  that  tho  design  criteria  required. 

The  major  difficulty  in  the  analysts  was  to  find  reliable  forms  of 
the  derivatives  in  terms  of  the  planfonr.  parameters.  The  forms  that 
follow  are  based  on  the  assumption  that  tha  effects  of  the  wing  and  body 
can  be  separated  and  added.  This  assumption  remained  valid  throughout 
the  range  of  linear  aerodynamics  considered  In  the  analysis*  Throughout 
tho  analysis  where  the  derivatives  depend  on  parameters  that  are  defined 
only  graphically  In  the  references,  those  derivatives  wl II  be  expressed 
as  functional  relationships.  Where  static  margin  was  a parameter,  a value 
of  .25  was  used  for  the  derivative  value  shown, 

C. 

a.  This  derivative  represented  the  change  In  lift  with  changing  angle 
of  attach.  Tha  theoretical  value  of  2*  was  assumed  and  corrected  for  finite 
wing  effects  using  a correction  factor  presented  In  reference  6. 

Cl  * ?(V  *R>  > 

Q 

C - 4.25368 
Si 

<*Q 

q.  This  derivative  ropresonted  the  change  In  drag  with  changes  In  angle 
of  attack. 


Cn  - f<CL,  CL  »®,AR> 
a o 


.5091 


(Ref.  7) 


(12) 
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Sy| 

a . This  dor  I votive  roprosented  tho  change  In  the  momont  coefficient 
with  changes  In  angle  of  attack.  It  represented  the  primary  factor  In 
static  stability  and  as  such  had  to  be  negative  to  ensure  static  stability. 


* *kn  Cl 
a 0 


(Ref.  4)  (13) 


a 


-1.05842 


-i  9m 

. ,&»  »hd  a.  These  derivatives  wore  assumed  to  be  negligible  In  a tail- 
less design  since  the  primary  contribution  to  the  magnitudes  of  the  der- 
ivatives was  a result  of  the  lag  In  downwesh  ranching  the  tall  (Ref.  4) 
and  thus  was  absent  In  this  configuration. 

C, 

q.  This  derivative  resulted  from  the  change  In  lift  due  to  pitching 
velocity. 


C - f(kn,  CL  ) 
q a 

C - 4.23368 


(Ref.  7)  (14) 


q.  This  derivative  resulted  from  tho  change  In  pltchlnq  moment  due  to 
changes  In  the  pitching  velocity.  It  represented  the  primary  damping 
force  for  short  period  oscillations. 


Cm  ■ f(A,  C , kn,  AR) 
q « 


(Ref.  8)  (15) 


X 


-.993186 


Tho  equilibrium  flight  condition  used  In  the  stability  analysis  was: 


1 


C « .86815 
Le 


CD  » .1200 

ue 


Y » -7.9* 
’e 
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The  characteristic  matrix,  eigenvalues,  and  eigenvectors  were  calculated 
by  a digital  computer  program  written  for  this  purpose  and  Argand  diagrams 
of  the  eigenvectors  are  presented  In  Appendix  A.  A root  locus  of  the 
eigenvalues  for  a range  of  static  margin  from  .05  to  .30  for  the  short 
period  is  presented  in  figure  5a  and  for  the  phugoid  in  figure  5b.  It  can 
be  seen  that  the  longitudinal  modes  were  both  stable  and  thus  represented 
satisfactory  solutions  to  longitudinal  stability. 


Fig.  5a  Root  Locus  Wing-Body  - Short  Period 
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Utilizing  the  same  equilibrium  conditions  as  In  the  longitudinal  case,  the 
lateral  derivatives  were  calculated  In  the  following  manner, 

Cy 

6.  This  derivative  represented  the  change  In  side  force  due  to  sideslip. 

The  major  contributions  were  from  body  effects  with  small  effects  from 

the  wing.  The  best  representation  for  this  derivative  was  created  by 

summing  dihedral  and  body  effects  from  reference  8 with  the  wing  effects 

presented  in  reference  7, 

C B f(r,SQi  S,  Ci  ,A  ,AR) 

Yfl  0 e (17) 

C = .324182 
v8 

Although  typical  values  of  this  derivative  are  negative,  these  values 
arise  from  tall  effects  which  overcome  the  destabilizing  effects  of  the 
body  and  are  not  present  In  this  configuration. 


C 


8.  This  derivative  was  of  paramount  concern  In  this  design  since  It 
was  the  major  determinent  of  directional  stability;  one  of  the  critical 
design  points.  The  major  Influences  In  this  derivative  were  from  wing 
sweep,  dihedral,  and  aspect  ratio.  The  form  used  for  this  derivative  was 
a summation  of  body  and  dihedral  effects  from  reference  8 with  the  wing 
effects  from  reference  7. 


C.  = f CA,  r ,AR,  b,  d,  CL  , zw)  (l8) 

8 e 

C.  = -.12885 

■ A 
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C 

ng.  This  derivative  provided  the  yaw  stiffness  to  the  vehicle  and  a 
such  had  to  be  positive  or  the  vehicle  would  have  been  statically  un- 
stable in  yaw.  Normally  the  effects  of  the  body  would  have  outweighed 
the  effects  of  the  wing  because  of  a large  body  area  ahead  of  the  center 
of  gravity;  however,  In  this  case  the  wlnq  effects  outweighed  the  body 
effects  and  caused  the  value  to  favor  static  stability.  Contributions 
from  the  body  were  evaluated  using  the  methods  of  reference  8 while  wing 
effects  were  calculated  from  methods  of  reference  7. 


» f(S, 


^B»  AtAR,  kn) 


-3 

8.8073x10 


(19) 


yD . This  derivative  represented  the  change  in  the  si  deforce  with  the 
change  in  wing-tip  helix  angle  (Ref.  3).  The  wing  contributed  the  major 
effects  to  this  derivative  and  since  experimental  evidence  indicated  that 
the  body  effects  were  negligible  (Ref.  3),  only  wing  effects  were  cal- 
culated. 


C - f(CL,  AR,  A) 
VP 

C « .158657 

Yn 


(Ref.  7)  (20) 


np.  This  derivative  represented  the  yawing  moment  caused  by  roll  and 
caused  the  close  coupling  of  roll  and  yaw  (Ref.  3).  In  the  case  of 
this  configuration,  these  effects  were  more  closely  coupled  than  normal, 
which  made  any  lateral  Instabilities  difficult  to  analyze  and  correct. 
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In  this  case  experimental  evidence  allowed  the  effects  of  the  body  to  be 
neglected  (Ref.  3). 

Cn  ■ f(AR,  A,  kn,  CL>  (Ref.  7)  (2|) 

P 

C - -.14163 
n 

p 

It  should  be  notsd  that  C varies  with  C.  as  shown  in  reference  3 and 

np  l 

therefore  could  have  been  of  either  sign.  The  negative  sign  in  the  pre- 
sent case  was  destabilizing  and  no  method  of  changing  Its  sign  could  be 
found  without  flying  at  a CL  that  was  unacceptable  from  a performanc 
consideration. 


v 

r.  This  derivative  represented  the  change  in  si deforce  with  variations 
In  yawing  velocity. 


C * , A,  AR,  kn)  (Ref.  7)  (22) 

yr  L 
C - -.03324 

' M* 


r.  This  derivative  represented  the  effects  of  asymmetrical  lift  and  drag 
distribution  over  the  wing  during  yaw  (Ref.  3).  In  the  present  case  the 
effects  of  the  body  were  possibly  significant  (Ref.  3),  but  no  known 
method  of  calculating  these  effects  was  known  (Ref.  3).  The  value  of  this 
derivative  was  similar  to  valuos  in  reference  9 so  the  value  was  assumed 
to  be  approximately  correct. 


C * f(C,  , AR,  k , Cn  , A)  (Ref.  7)  (23) 

n...  l n ue 

Cn*  = -.030108 
r 
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vi 

p.  Thl*  derivative  represented  the  change  In  rolling  moment  coefficient 
with  the  change  In  wlng-tlp  helix  angle  and  In  this  cm«  experimental 
evidence  Indicated  that  body  effects  were  negligible  (Ref.  3). 


C.  ■ -.405159 
P 


(Ref.  7)  (24) 


Jr.  This  derivative  represented  the  rolling  moment  due  to  yaw  rate. 
Fuselage  effects  could  also  be  neglected  in  this  case  (Ref.  3) 


C,  ■ .23995 
r 


(Ref.  7)  (25) 


The  equilibrium  conditions  for  the  lateral  analysis  were  the  same  as 
those  used  In  the  longitudinal  analysis.  A root  locus  for  the  dutch 
roll  for  a static  margin  range  from  ,05  to  .30  is  shown  In  figure  5c. 
Argand  diagrams  for  the  static  margin  range  are  shown  In  Appendix  B. 

As  shown  In  the  root  locus  below,  the  vehicle  was  found  unstable  in 
dutch  roll.  No  method  was  found  to  stabilize  the  analytical  model  of  the 
vehicle  from  changes  In  the  vehicle  parameters,  such  as  wing  configur- 
ation, mass,  moments  of  Inertia,  or  bony  shape. 
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Wind  Tunnel  Tests 

Wind  tunnel  tests  were  performed  In  an  attempt  to  validate  the 
analytical  model,  find  possible  corrective  measures  for  the  apparent 
lateral  Instability,  and  to  provide  Information  In  areas  not  covered  by 
the  analytical  model.  The  tests  were  performed  In  the  AFFDL  one-meter 
low  speed  wind  tunnel;  a closed  loop,  constant  speed  tunnel.  All  data 
runs  were  conducted  at  the  tunnel  airflow  maxlnum  speed  of  75fps., 
which  closely  simulated  the  analytical  model's  84  fps.  Runs  were  made  to 
measure  both  lift  and  drag,  with  moment  data  unavailable  due  to  equip- 
ment malfunction.  Curves  showing  the  data  are  presented  in  figure  6. 
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Fig.  6 Wing-Body  Lift  and  Drag  Curves 


As  seen  In  figure  6,  both  CL  and  Cq  are  within  20J»  of  that  predicted  by 
the  theoretical  model.  Three  further  studies  were  performed  to  Investi- 
gate model  performance  and  to  partially  validate  the  theoretical  model. 

The  first  of  these  runs  was  performed  after  tufting  the  model  to 
Investigate  tho  airflow  around  the  model.  This  test  was  performed  at 
.agles  of  attack  of  I0#,  15*,  and  20*  to  validate  the  analytical  model 
In  flight  conditions  around  the  equilibrium  position.  Very  smooth  air- 
flow was  observed  over  the  surface  of  the  entire  vehicle.  In  particular 
around  the  wing  body  juncture.  This  area  of  smooth  f*ow  helped  to  vali- 
date the  assumption  that  body  and  wing  effects  could  be  calculated 
separately  and  added  together. 

The  second  run  was  performed  in  a one-degree  of  freedom  mount  (Fig. 7) 
to  test  longitudinal  stability. 
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Fig*  7 One  Degree  of  Freedom  Long  I tud Inn I Mount 

This  run  enebled  both  the  trim  angle  of  attack  to  he  found  for  this  dynam- 
ic pressure  and  the  reflex  angle  chosen  for  the  camber  line  to  be  tested* 
As  the  airflow  was  started,  the  model  trimmed  at  an  angle  of  attack  of  12* 
and  was  vory  stable  In  pitch.  Motions  induced  by  Interruptions  In  the  air- 
flow were  damped  out  rapidly  as  predicted  by  tho  short  period  eigenvalue 
analysis*  In  this  case  both  the  reflex  angle  predicted  and  the  short 
period  response  were  validated  by  the  experimental  evidence. 

The  third  test  was  conducted  In  a similar  manner  to  that  of  the  short 
period  test  with  the  exception  that  the  model  was  mounted  In  a fitting 
that  allowed  freedom  around  the  z-axis  (Fig.  8).  Although  friction  In  the 
mounting  bearing  could  have  effected  the  damping,  it  could  not  have  pro- 
vided a restoring  force  so  this  run  demonstrated  the  directional  stabil- 
ity of  the  vehicle,  since  the  model  trimmed  straight  ahead  and  any  per- 
turbations wore  damped  out  rapidly.  None  of  the  closely  coupled  modes 
of  the  analytical  model  were  ablo  to  bo  completely  validated  due  to  a 
lack  of  sufficient  degrees  of  freedom  In  the  tunnel  mounts. 

The  major  conclusions  provided  by  the  wind  tunnel  tests  were  the 


24 


GAE/MC/'MD-I 


lift  and  drag  data,  and  partial  validation  of  tha  analytical  modal.  In 
particular,  tha  longitudinal  short  parlod  mode  was  validated  as  was  tha 
airfoil  design,  No  fasts  could  ba  davlsad  with  tha  aval  labia  equipment 
to  tast  tha  unstabla  lataral  moda,  so  no  conclusions  could  ba  mada  con* 
earning  total  vahlcla  stability  or  mathods  to  corract  tha  lataral  Insta- 
bill ty  pradlctad  by  tha  analytical  modal. 


Airflow 

Into 

Pepar 


Fig.  8 Ona  Degree  of  Freedom  Lateral  Mount 


Moment  o?  Inertia  Calculation 

The  following  method  of  moment  of  Inertia  calculation  was  discussed 
extensively  In  reference  10,  For  inertia  measurement  about  the  y-exis, 
a pivot  was  inserted  Into  the  model  fuselage  at  right  angles  to  the 
x-axIs  and  vertically  above  the  longitudinal  center  of  gravity  position. 
Tha  modal  was  then  suspended  from  this  olvot  at  a distance  I from  the 
canter  of  n-evlty  so  the  model  oscillated  about  an  axis  parallel  to  the 
y-axls.  The  moment  of  Inertia  about  the  y-axls  was  then: 

| . Twl  - wl3  (26) 

vv  TO  g 
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•yy  was  also  a principal  moment  nun  1o  tho  x-z  plana  of  svnmetry.  Tha 
partod  was  determined  to  an  accuracy  of  , It  by  timing  50  oscillations 
with  a stopwatch. 

To  determine  the  principal  axes  orientation,  the  model  was  set  at  a 
series  of  small  pitch  angles  to  the  horizontal  and  the  rolling  moment  of 
Inertia  was  determined  at  oach  attitude  by  the  above  method.  The  minimum 
rolling  moment  of  Inertia  obtained  In  this  manner  was  the  principal 
rolling  moment  of  Inertia  and  the  angle  with  the  horizontal  gave  the 
orientation  of  the  principal  axes.  The  rosul+s  of  this  calculation  Indi- 
cated that  the  y-z  plane  was  also  mass  symmetric,  therefore  all  three 
axes  were  principal  axes. 

The  yaw l no  moment  of  Inortfa  was  determined  by  swinging  the  model  as 
a torsional  pendulum,  using  i bt-filar  suspension.  The  two  wires  were 
attached  to  the  x-exls  on  egual  distance  fore  and  aft  of  the  center  of 
gravity.  Then  the  yawing  moment  was: 


1 


zz 


wT*d 

TOT 


C7) 


These  Inertias  were  corrected  for  confer  of  gravity  location  for  the 
various  static  maroins  ar.d  non- dimensional !zed  using  the  non-d I mens lon- 
al  izatlon  from  reference  4 so  that  they  would  be  compatible  with  those 
eguatlons.  For  a static  margin  of  .25  tho  Input  Inertias  wero; 

I ' * 33.209 
xx 

iyy*  - 770.120 

? * - 35.11? 

zz 
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Might  Tost 

This  portion  of  the  analysis  was  performed  by  hand  launch  from  a 
ha 1 1 cop tar  moving  at  a forward  valoclty  of  60  knots  and  flying  at  an 
latitude  of  790  ft,  Thasa  tasts  showad  that  tha  vahleta  was  vary  un- 
stabla.  In  what  appaarad  to  ba  tha  dutch  roll  mode,  although  this  ob- 
sarvatlon  was  vary  subjactiva.  Slnca  thasa  tasts  conflrmad  that  tha 
vahlcla  was  unstabla  as  pradlctad  by  tha  analytical  modal,  soma  typa  of 
vahlcla  modification  was  raqulrod  to  achlava  stability, 

Naw  Conf I curat Ions 

Thraa  typos  of  modifications  wara  triad  In  an  attampt  to  stabilize 
tha  vahlcla.  Tha  first  of  thasa  was  tha  addition  of  a vortical  stabll- 
I rar  to  tha  raar  of  tha  modal.  This  configuration  was  tunnal  tastad  and 
It  appeared  to  damp  out  tha  parturbatlons  In  sldaforce  mora  rapidly  than 
tha  basic  vahlcla.  Whan  this  vahlcla  was  flight  tastad  vary  little  Im- 
provement was  noted,  so  without  a longer  tall  moment  arm  this  modifica- 
tion was  Insufficient,  Tha  next  two  attempts  wore  similar  In  Intent. 

In  an  attampt  to  prevent  tha  vehicle  from  entering  the  oscillations, 
two  type',  of  drag  devices  wara  triad.  Tha  first  was  a ribbon  parachute 
2"  wide  and  12"  long.  Tha  second  device  was  a conical  rigid  parachute 
2"  In  dlanatar  attached  by  a 12”  thread  to  tha  raar  of  the  vehicle  on  an 
axis  through  the  canter  of  gravity.  In  both  cases  the  additional  drag  was 
so  large  that  tho  weight  component  was  insufficient  to  provide  flying 
airspeed. 

Conclusions 

After  the  last  tests  of  the  vehicle.  It  was  decided  that  this 
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conf Iquratlon  would  not  be  further  Invast  I gated.  It  was  unllkaly  that  a 
•tibia  form  of  this  configuration  could  ba  developed,  and  avar  If  marqlnal 
stability  wars  achieved,  It  would  not  ba  as  satisfactory  a configuration 
as  tha  conventional  vehicle. 
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IV.  THE  WHKMiQPY-TAt L DESIGN 

As  major  stability  problems  began  to  appear  In  tha  wing-body  analy- 
sis, mora  emphasis  was  placad  on  tha  w I no-body- ta 1 1 daslgn  study.  Tha 
first  step  In  tha  das  I on  process  was  to  define  In  exact  terms  tha  plan- 
form  parameters  for  the  first  Iteration  of  tha  wing- body- to 1 1 analysis. 

Tal I Dos  Ion 

Tha  best  location  from  stability  considerations  of  the  tail-boom 
body  Juncture  was  high  on  tho  body  (Ref.  II).  This  study  of  a similar 
configuration  found  that  tha  vehicle  had  greater  lateral  stability  with 
tha  tail-boom  located  high  on  the  body,pr Imarl ly  due  to  the  elimination 
of  body-tall  Interference  at  some  high  angles  of  attack.  Therefore,  the 
tel  I -boom  was  mounted  Just  below  tha  wing  trailing  edge  as  shown  In  fig- 
ure 2.  It  was  seen  In  the  wind  tunnel  tests  of  this  configuration  that 
this  location  did  eliminate  body-tall  Interference  as  desired. 

The  specific  details  of  the  tall  design  were  arrived  at  from  an 
evaluation  of  the  packaging  problem.  The  first  tal I -boom  folding  mechan- 
ism considered  was  a pivot  at  the  body-boom  juncture  with  the  boon  being 
constructed  from  a material  with  sufficient  flextbllltv  to  allow  It  to  be 
wrapped  around  the  circumference  of  the  body.  The  only  practical  ver- 
tical tall  planform  useoblo  with  this  tal I -boom  design  was  a twin  vertical 
tall  with  the  vertical  stabilizers  located  at  the  extreme  onds  of  tho  hor- 
izontal stabilizer.  To  provide  the  maximum  available  horizontal  stabil- 
izer span,  tho  span  was  chosen  so  that  tho  vertical  stabilizers  would 
fit  outside  of  the  wings  whon  folded. 

Tall-bocn  lonqth  was  choson  as  a compromise  between  boom  flexibility 
and  tall  sensitivity.  Since  the  toil-boom  was  to  be  flexible,  the  shorter 
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the  tail-boom  within  the  Units  of  stability,  the  less  sensitive  the 
vehicle  would  be  to  tail  datamation.  For  this  reason,  a tail-boom  length 
of  9”  was  chosen  as  the  starting  point  In  the  design  process.  Both  hori- 
zontal stabilizer  chord  and  vertical  tail  area  were  chosen  to  provide  a 
tal!  volume  of  slightly  over  .5  since  the  tall  volume  of  the  configuration 
In  -eference  II  was  approximately  thjs  value. 

Stabi i i ty  Ana  lysis 

The  same  assumptions  made  in  tne  wing-body  analysis  were  necessary 
in  developing  this  analytic  model,  except  where  less  stringent  assumptions 
were  allowable  as  noted  in  the  analysis.  One  further  assumption  was  made 
concerning  flow  at  the  tail.  It  was  assumed  that  the  tall  w^s  located  in 
an  area  of  undisturbed  *low  except  for  the  downwash  effects  on  the  tail 
angle  of  attack.  The  evidence  of  reference  II  Indicated  that  this  was  a 
valid  assumption.  The  same  equations  of  motion  and  the  same  resulting 
characteristic  matrix  as  used  in  the  w inn-body  analysis  were  used  in  this 
analysis  and  therefore  will  not  be  repeated. 

Longitudinal  Analysis.  The  equilibrium  conditions  used  in  this  analysis 
were  similar  to  those  used  in  the  wing-body  analysis  and  are  not  repeated. 
Where  the  stability  derivatives  are  a function  of  static  margin,  a static 
margin  of  .25  was  chosen  for  demonstration  purposes.  All  of  the  following 
derivatives  and  those  In  the  lateral  analysis  were  calculated  using  the 
methods  of  reference  8. 

CL  = f (AR,  SH,  C,  , C.  ) (28) 

° aw  °H 
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Cl  - 4.524 

°Theory 


4.756 


“wind  tunnel 


The  velue  of  2 n Mas  assumed  as  the  horizontal  stabilizer  C.  because  of 

'a 

the  endplate  effects  of  the  vertical  stabilizer.  As  seen  from  the  wind 
tunnel  test  data  the  resulting  value  was  quite  good  using  this  assumption. 


CD  * f(CL»  CL 
o 

CQ  « .64848 

aTheory 


u,lnd  tunnel 


.1203 


AF) 


(29) 


The  predicted  value  was  not  very  good  In  this  case.  Since  the  forms  of 
the  derivatives  used  In  this  analysis  were  empirical  curve  fits  to  conven- 
tional data  sets.  It  Is  possible  that  the  lack  of  a conventional  fuselage 
resulted  in  an  artificially  high  value  being  predicted. 


S)  = * kn^  (30) 

a a 

Cu  = -1.1890 


C = f(C|  , C,  , kn,  VH>  (31) 

«q  aw  aR 

C,  = 9.8775 

q 
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V 

—3.* 


Si  “ f<Cl  * Cl  * kn»  VH*  AR> 

q «w  aH  (32) 


a.  - -19.157 

^q 


a. 

tall 


The  value  for  this  derivative  could  only  be  calculated  for  the 
effects  and  was  thus  multiplied  by  a factor  of  l.l  (Ref.  4). 


V 


-4 


"C,  ,VU>  fe 


H*  da 


(33) 


2.5222 


JM. 

a. 


s ■ V fr 

C - -7.115 

Ma 


(34) 


The  value  of  the  pitching  moment  was  calculated  by  the  same  methods  used 
for  the  wing  body. 

t * 2499.4 

y 

The  resulting  eigenvalues  and  eigenvectors  for  a static  margin  of  .25  : 


Short  Period: 

-6,7864x1 0"3  * 2. 1 I98xl0"2 i 

C=  .3049 

A 

V: 

2.6008x10"?  * 1 . 5650x1 0”3 i 
4.5040x10  * 5.5652x10"' i 

a: 

A 

q: 

1. 021 5x1 0"2*  l.l 705x10”? } 

0: 

3.6909x10"'  * 5.9744x10"' i 

Phugoid: 

-5,5226x10"® * 7.7737xl0"4i 

C = .0709 

A 

V: 

4.5196x10” ' * 3.7087x10"' i 

a: 

-8.6222x10”^  * 5. 3903x1 0“3 i 

5.0104x10":  * 3.8555x10"? i 

0: 

-5.3904x10  *6.0623x10"' i 
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The  root  locus  for  a static  margin  range  of  ,05  to  ,30  for  the  short 
period  is  shown  in  figure  9a  and  for  the  phugoid  In  figure  9b.  Argand 
diagrams  for  the  static  margin  range  are  presented  in  Appendix  C.  As 
shown,  the  vehicle  was  stable  In  both  longitudinal  modes;  however,  the 
phugoid  was  weakly  damped  so  further  studies  were  made  of  the  sensitivity 
of  the  analytical  model  to  changes  In  the  derivative  values. 


Fig.  9a  Root  Locus  for  Wlng-Body-Tal I - Short  Period 
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Fig.  9b  Root  Locus  for  Wlng-Body-Tal  I - Phugold 


Longitudinal  Sensitivity  Analysis,  Two  results  were  obtained  from  this 
analysis:  I)  the  sensitivity  of  the  model  to  errors  In  the  derivatives 
was  evaluated,  and  2)  possible  methods  to  increase  phugold  damping  were 
found.  The  changes  In  each  derivative  and  the  resulting  eigenvalues  are 
shown  In  Appendix  0. 

It  was  seen  from  this  analysis  that  the  model  was  relatively  insen- 
sitive to  derivative  errors  since  only  small  variations  in  response  oc- 
curred and  the  model  remained  stable  for  all  variations.  The  only  deriv- 
ative that  made  noticeable  changes  in  the  phugold  was  (other  than 

q 

which  could  not  do  easily  changed)  and  even  this  derivative  made  only 
minor  changes  in  the  phuc  Id.  However,  as  a result  of  this  possible  im- 
provement in  response,  two  different  attempts  were  made  to  Increase 
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phugold  damping.  The  first  change  was  an  increase  In  the  horizontal 
stabilizer  chord  to  3M.  The  second  change  was  to  Increase  tail  volume  to 
.8,  which  although  similar  to  changing  the  chord,  could  be  arrived  at  by 
other  means  Including  lengthening  the  tail-boom,  or  Increasing  the  hor- 
izontal K.i  span.  Neither  of  those  two  changes  made  a noticeable  change 
In  the  phugold  response. 

One  completely  different  analysis  was  made  by  assuming  that  the 
phugold  response  was  simply  an  energy  transfer  between  the  potential  and 
kinetic  energies  (Ref.  4).  Using  this  representation  of  the  phugold.  It 
was  possible  to  see  one  further  change  that  would  definitely  Increase  the 
phugold  damping.  The  change  was  to  Increase  the  drag  of  the  vehicle, 
which  although  undesirable  from  a performance  standpoint,  could  easily  be 
adopted  If  necessary.  The  results  of  increasing  the  drag  coefficient  to 
,12  and  reducing  the  lift  coefficient  to  .7  to  Increase  the  flight  path 
angle  so  that  the  weight  component  would  balance  tie  Increase  In  drag  were: 


Short  Period: 

-6. 7859x1  O’3 

2. 1 195x10“^  I 

C » .3049 

0: 

2. 5799x1  O’3 
4.4371x10 

I.0985xl0"3l 
5,6175x10“! 1 

a; 

1.0353x10““ 

3.5384x10 

1.1583x10*71 
6.0175x10'*  1 

0: 

Phugold: 

-I.l523xl0"4 

6. 2474x1  O'4 1 

C = .18064 

0: 

5.7758x10'* 

8.8540x10**^  i 

a: 

-8.2241x10  j 

2,6877x10.  i 
2.4096x1 0~?  1 
7.9044x10  1 

« 

q: 

5.1707x10  , 
-1,8358x10'* 

e: 

As  shown,  this  change  resulted  in  a major  change  in  the  phugold  response; 
however,  the  glide  ratio  was  reduced  by  50J.  to  5.8.  This  glide  ratio  was 
acceptable  to  meet  the  performance  regu 1 remonts  If  that  large  of  a trade- 
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off  In  performance  was  necessary  to  obtain  the  desired  response.  As  a 
result  of  this  study*  It  was  decided  that  the  oriqinal  phugold  response 
would  be  retained  unless  flight  tests  showed  that  It  was  unacceptable. 

Aeroelastlc  Analysis.  The  influence  of  tail-boor*  flexibility  was  consi- 
dered to  have  a major  impact  on  vehicle  stability  because  of  the  proposed 
tall  packaging  concept.  The  following  analysis  was  conducted  under  the 
assumption  that  the  tall  flexibility  could  be  represented  by  assuming 
that  a torsional  spring  located  at  the  pivot  provided  all  of  the  flex- 
ibility to  the  boor**  while  the  rest  of  the  boom  was  assumed  to  be  rigid. 
This  assumption  was  valid  because  the  only  significant  change  occurring 
because  of  tall  flexibility  was  a change  In  the  angle  of  att ick  of  the 
tall  relative  to  the  vehicle*  and  no  significant  effocts  were  caused  by 
the  boom  itself.  The  following  analysis  follows  that  of  reference  4 
closely.  Assume  that  Vo,  * -kLt  where  k is  the  flexibility  of  the  tor- 
sional spring  located  at  tho  pivot.  The  figure  below  Illustrates  the 
assumed  model. 


Fig,  10  Aeroelastlc  Model 
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let: 


*t  * CL,f  (°W0  * C * V 


(35) 


then: 


* at  rlald 

■f  I + kt.’  l.5eV^J 

at  rigid 


(/,WQ  * c * 


(36) 


and  the  tall  effectiveness  has  baan  raducad  by  a factor  of: 


t r v? 

marL ctpTO 

at  rigid 

After  this  factor  was  substituted  Into  the  appropriate  places  In  the 
stability  analysis,  the  following  reductions  In  static  margin  wore  found 
for  the  flexibility  values  shown  at  the  basic  static  margin  of  .25: 

kbaslc  " (flight  test  model  flexibility) 

Akn  ■ .06 

k-  .08818  (1.5  kbaslc> 

Akn  ■ .08 

k • .03937  (keas|c  / 1.5) 

&k  ■ .04 
n 

The  changes  in  the  static  margin  were  significant  for  the  case  where 
either  pay I os'  mqu»=  nets  required  a low  static  margin  to  meet  the 
packaging  requirement,  or  the  choice  of  a material  for  the  tail-boom 
had  more  flexibility  than  evaluated  In  the  above  analysis  thus  resulting 
In  a greater  decrease  static  margin  than  predicted.  The  solution  to 
this  problem  will  be  further  discussed  In  the  Iteration  phase. 
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Lateral  Stability  Analysis.  The  same  equilibrium  conditions  usod  In  tho 
longitudinal  analysis  were  usod  In  tho  following  analysis* 


where: 


cy  - f«s0.  s, 

5 


t 

-.80262 


(38) 


C,6  * f<CLWB*  AR*  ,v*  Cy0  * A'  0>  (39) 

C « -.12442 


C * ( ( S,  I , I , z , C i <*) 

"0  8 V V V0  (40) 

v 

C - .05487 
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■ f(zv. 


I v» 
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C « .02496 
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Ci  * ,(bH#  Zv'  SH* 
P 

C,  - -.42050 
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Cn  * MCI  * CL»  *V  V CV(I  * a>  U3> 

C • -.12461 
n 


C . «IV,  *v.  C . •> 

Tr 

C - .44660 
“r 


(44) 


Cl  * f(CL»  ‘v*  V Cy  * 0,)  <45) 

f By 

C » .22013 
*r 


f<cL.  CQ  . 

o 


- -.26046 


(46) 


The  resulting  eigenvalues  tor  this  analysis  were: 

Dutch  Roll  Rolling  Convergence  Spiral 

4.0907xl0’4  * 0. 1 1 72x1 o”3 I -4.4639XI0-3  -2.7798x10’ 

The  resulting  root  locus  tor  the  dutch  roll  mode  is  shown  In  tlgure  II. 
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Arqand  diagrams  ♦or  the  static  margin  range  are  presented  In  Appendix  E. 
As  shown  In  the  response,  the  dutch  roll  was  unstable.  To  find  possible 
solutions  to  the  Instability  other  than  an  excessively  large  static 
margin  which  would  have  required  the  addition  of  wlnq  sweep  and  to  eval- 
uate model  sensitivity  to  derivative  errors,  a sensitivity  analysis  of 
the  effects  of  the  variation  of  the  Individual  derivatives  on  the  lateral 
response  was  made. 

Lateral  Sensitivity  Analysis.  The  amount  of  derivative  variation  and  the 
resulting  eigenvalues  Is  shown  in  Appendix  F.  The  results  of  this  analy- 
sis showed  that  the  error  sensitivity  was  minimal  oxcept  in  the  case  of 

C,  and  that  this  derivative  was  the  only  derivative  within  the  variation 
*6 
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studied  that  caused  vehicle  stabilization. 

The  result  obtained  from  the  sensitivity  analysis  was  confirmed  by 
the  approximate  solution  to  the  dutch  roll  (Ref,  4) 


DR 


1/2 


approx 


"OR  . -4 

approx  ■ 4.1736x10 

A 

Rnp  «4 

analyt  • 4,0907x10 


This  approximation  also  Indicated  that  increases  ,n  cyg  • cn  » an<J  Cna 
would  result  In  a stable  dutch  roll  mode.  After  Investigating  the  func- 
tional dependence  of  each  of  the  concerned  derivatives.  It  was  determined 
that  a change  In  vertical  stabilizer  chord  would  cause  the  desired  varia- 
tion In  all  of  the  necessary  derivatives.  Therefore,  in  Increase  In  chord 
to  3"  was  evaluated.  The  resulting  eigenvalues  and  eigenvectors  were: 


Dutch  Roll  Rol I Ing  Converoenc*  Spiral 


-2.6063x1  O’4 

8.91 79xl0"3 1 

-3.8532x1 0*3 

-2,6755x10*“ 

3.5547x10“' 

4.5351x10*' 1 

-4.64l3xl0~* 

1.1013x10* 

1.9794x10”) 
-1.9655x10  : 
6,5921x1 0”1 

3.0243x10“; 1 

-2.0062xl0”f 

-1. 1998x10. 

1.8271x10"“  1 
4,71 1 1x10"  1 

2.0508x10  ' 
9.9872x10"' 

2.9588x10 

9.9993x10" 

The  resulting  root  locus  for  the  entire  static  margin  range  is  shown  In 
figure  12. 
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Fig.  12  Root  Locus  for  Wlng-Qody-Tal  I - Dutch  Roll  2 


The  resulting  response  was  stable  as  desired  and,  although  the  dutch 
roll  was  of  long  period.  It  was  satisfactory  to  neot  the  design  re- 
quirements. The  response  for  the  static  margin  range  Is  presented  In 
Appondlx  G. 

The  effects  of  altitude  which  were  Important  In  this  design  were 
discussed  at  length  In  reference  4 and  since  they  mainly  enter  through 
the  non-dimensional Izotion  of  the  terms,  they  were  not  studied  explicitly. 
As  discussed  In  reference  4,  the  effocts  of  altitude  variation  were  main- 
ly to  reduce  damping  as  altitude  was  Increased  and  to  Increase  the  period 
of  both  the  phugold  and  the  dutch  roll  modes.  Although  those  effects  were 
undesirable,  the  only  correction  for  these  effects  was  to  make  the  final- 
ized vehicle  as  heavily  damped  as  possible  in  these  two  modes  without 
sacrificing  stability  in  the  remaining  modos. 
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Wind  Tunnnl_  l£*li 

The  MM  set  of  wind  tunnel  tests  were  run  for  the  wing-body- tel  I 
configuration  at  were  run  for  the  wing-body  configuration.  Both  lift  end 
drag  were  measured,  and  the  other  tests  to  evaluate  the  analytical  model 
end  the  assumptions  necessary  in  Its  development  were  run.  The  lift  end 
drag  curves  that  resulted  are  shown  In  figure  13. 


The  airflow  tests  In  this  case  were  accomplished  through  the  use  of  a 
bubble  generator  which  dispersed  noutrnlly  bouyant  jubbles  Into  the  air** 
stream.  This  method  allowed  better  evaluation  of  the  airflow  around  the 
tall  than  was  possible  with  tufting.  The  flow  around  the  body  was  expected 
to  bo  very  similar  to  that  around  the  body  of  the  wing-body  vehicle  since 
the  bod  I os  were  Identical,  and  thus  would  not  nave  to  be  studied  as 
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carefully  as  tufting  would  have  allowed.  The  flow  around  the  tall  was  of 
major  Interest  since  any  tall  blockage  within  the  normal  angle  of  attack 
range  would  have  Introduced  some  type  of  Instability  Into  the  model  dy- 
namics and  would  also  have  Invalidated  the  analytical  model  In  the  affect- 
ed range.  The  resulting  streamlines  showed  no  tall  blockage  at  any  angle 
of  attack  from  -10°  to  +20°. 

In  this  case  as  In  the  tailless  case  only  the  short  period  mode  and 
the  directional  stability  could  be  evaluated  using  the  one  degree  of  free- 
dom mounts.  In  both  cases,  the  modes  were  very  stable  and  well  damped. 

FI Ight  Tests 

The  flights,  conducted  In  the  same  manner  as  the  wing-tody  tests, 
resulted  In  demonstrating  that  the  model  was  stable  In  ^11  modes.  The 
major  problem  that  appeared  could  have  been  either  slow  damping  of  the 
launch  perturbations  from  the  helicopter  rotor  wash,  or  a wealtly  damped 
dutch  roll  mode  with  large  amplitude  convergent  oscillations  that  resem- 
bled an  exaggerated  dutch  roll  oscillation.  This  problem  could  not  be 
further  evaluated  In  the  flight  test  program  since  the  vehicles  were  not 
of  sufficient  size  to  see  at  an  altitude  that  would  allow  flight  dura- 
tions of  sufficient  length  to  see  these  perturbations  damp  out.  Although 
this  design  did  meet  all  of  the  design  criteria,  certain  problems  still 
existed  that  a further  iteration  of  the  design  process  would  solve. 

Design  Iteration 

Three  major  problems  needed  to  be  solved  if  the  design  was  to  be 
very  successful:  I)  the  tail  packaging  problem  needed  a batter  solution 
so  that  the  trade-off  between  packaging  and  stability  would  not  be  as 

great,  2)  the  tail  flexibility  problem  had  to  be  eliminated  so  that  less 
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limitation  on  payioad  center  of  gravity  location  was  necessary,  and 
3)  the  oscillations  In  the  dutch  roll  mode  as  predicted  by  the  analytical 
model  and  seen  in  the  fliqht  tests  needed  to  be  reduced.  The  first  two 
problems  could  be  reduced  or  eliminated  by  a different  toll  foldlnq  mech- 
anism, and  the  dutch  roll  mode  could  bo  Increased  by  an  increase  In  the 
vertical  tall  effectiveness  as  shown  In  the  sensitivity  analysis.  The 
tall  modification  presented  In  flqure  14  was  developed  to  solve  all  of 
the  problems. 


The  resultinn  tail-boom  was  designed  so  that  the  vertical  surface  folded 
flat  against  the  horizontal  surface  and  the  flat  surface  formed  then  wrap- 
ped around  the  circumference  of  the  body.  ..hen  unfolded,  the  vertical  and 
horizontal  surfaces  formed  a symmetrical  cross  that  was  rigid  In  the  hor- 
izontal and  vertical  planes,  thus  eliminating  the  flexibility  problem. 
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Although  torsional  rigidity  was  reduced  with  this  torn,  the  elastic  axis 
of  the  tail -boom  extended  through  the  axis  of  symmetry  of  both  the  hor- 
izontal and  vertical  surfaces  thus  creating  a symmetrical  load  distribu- 
tion In  symmetrical  flight.  Asymmetrical  loadings  could  not  be  analyzed 
using  the  available  model  and  therefore' pending  wind  tunnel  analysis  were 
assumed  negligible.  The  response  of  the  resulting  modification  using  the 
dimensions  Illustrated  In  figure  14  was: 


Long i tud i na I : 


Short  Period 

•3.9622x1  o’3  ± 2.0363xlO"Z  I 

■2.0747x1 0"4  ± 3. 5663?: I 0"3  I 
2.8357x10"'  * 6.5245x10“*  I 
*1.3605x10 "f  * 5.2340xl0"3  l 
3.7237x10"  5.9558x10“*  I 


Phugold 


-5.5374x10 

2. 5306x1  o' 
-3.6446x1 0_ 
3.0338x10' 
-7.0412x10' 


7.8522xl0_4i 

5.2263xlo“S 
4,3088xl0”ji 
5. 6603x10  , I 
3.2994x10  I 


Latera I : 


Dutch  Roll  Rolling  Mode 

■I.l909xl0“3  i.2965xl0“a  I -2.2526xlo"2 


1 . 1 339x1 0"^ 
-4.3672xl0“3 
4.8779x10"-; 
-3.l064xl0"2 


7.1615x10  ' l 
6.201 0xl0"3  I 
3.6520x1 0“3  i 
6.9405x10"'  i 


2.0282x10' 

-1.1620x10' 

2.9664x10' 

9.9993x10' 


Spiral 

-2.8964xl0"4 

-6.6366x1  O’*3 
I.3l68xl0“3 
-2.9440xl0“3 
-9.9997x10“' 


In  this  analysis  the  inertial  properties  were  assumed  duplicated  in  the 
modified  vehicle  due  to  lack  of  sensitivity  to  Inertia  variation  demon- 
strated in  the  sensitivity  analysis.  As  seen  in  the  root  locus  for  these 
modes  the  response  is  significantly  improved  over  the  baseline  vehicle. 
Figure  15  illustrates  the  root  locus  for  each  oscillatory  mode. 
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Fig.  15c  Root  locus  for  Cruciform  Tall  - Dutch  Roll 


Wind  Tunnel  Tests 

All  of  the  wind  tunnel  tests  run  on  the  other  two  configurations 
were  repeated  for  this  modified  version.  The  only  significant  differ- 
ence from  the  results  already  discussed  for  the  baseline  tallied  vehicle 
were  the  lift  and  drag  curves  presented  in  figure  16.  It  should  be  noted 
that  the  glide  ratio  improved  greatly  over  the  baseline  vehicle.  This 
Improvement  was  due  to  the  reduction  In  base  drag  which  was  a result  of 
the  new  tail-boom  aerodynamical ly  fairing  in  the  body  so  that  aerodynam- 
ical ly  It  seemed  longer  than  It  was  really. 
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Fig.  16  Cruciform  Tall  Llfl  and  Drag  Curves 


FI inht  Test  Analysis 

No  aircraft  support  was  available  to  provide  an  airborne  platform 
to  launch  the  vehicle  from,  so  this  set  of  flight  tests  had  to  be  hand- 
launched  from  tho  top  of  a building.  Since  there  was  no  similarity  In 
the  post-launch  disturbances  between  this  vehicle  and  the  previous  vehi- 
cles, direct  comparison  was  impossible.  Even  so,  It  was  easily  observed 
that  this  vehicle  was  more  stable  than  the  other  vehicles  tested  and  thus 
the  design  was  partial Iv  proven  by  this  sot  of  tests. 

Iteration  Conclusions 

As  a result  of  this  limited  set  of  tests  and  the  other  data  on  the 
configuration,  it  was  decided  lhat  this  vehicle  was  a satisfactory  design 
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to  submit  to  the  Avionics  laboratory.  The  next  step  In  the  analytic 
procedure  was  to  evaluate  the  final  vehicle’s  performance. 

Performance  Analysis 

In  order  to  evaluate  the  performance  of  the  vehicle  without 
arbitrarily  fitting  the  drag  to  a parabolic  polar,  the  method  of  ref- 
erence 12  was  used  In  the  performance  analysis.  It  was  shown  In  this 
reference  that  It  was  possible  to  evaluate  the  coefficient  of  lift 
and  thus  the  angle  of  attack  for  both  maximum  endurance  and  maximum 
range  from  a graph  of  CL  vs.  CD  In  the  first  instance  and  C^2  vs.  CD 
In  the  second  instance.  The  two  curves  are  shown  In  figure  17a  and  1 7b. 
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Fig.  1 7b  Maxlnun  Endurance  Calculation 

As  presented  In  the  reference,  the  point  of  maximum  range  occurs  at 
the  point  on  the  CL  vs.  Cq  curve  where  a line  drawn  from  the  origin 
Is  tangent  to  the  curve.  In  a similar  fashion,  the  point  of  maximum 
endurance  Is  found  on  the  second  curve.  The  performance  points  are: 

Maximum  Range:  « .8500 

« 8.5® 

Maximum  Endurance:  * 1.0000 

« 10.5° 

One  further  factor  had  to  be  taken  Into  account  before  the  trim 
angle  of  attack  could  bo  set.  At  speeds  below  the  speed  for  L/D 
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an  aircraft  Is  unstable  and  any  disturbance  will  cause  the  vehicle  to 
decrease  speed  and  to  eventually  stall*  therefore  the  trim  ancle  of 
attack  had  to  be  a compromise  between  speed  stability  and  maximum  ranee. 
Since  the  exact  function  of  the  payload  will  probably  vary  during  any 
production  life  of  the  vehicle*  a trim  ancle  of  7*  was  chosen  to  repre- 
sent e solution  close  to  both  the  optimum  angles  for  best  range  and 
best  endurance  and  was  still  a speed  stable  angle  of  attack. 
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V.  CONCLUSIONS  AND  RECOWENDATIONS 

Conclusion! 

The  major  conclusion  arrived  at  from  the  analysis  was  that  the 
cruciform  toll  design  was  the  best  design  to  submit  to  the  Air  Force 
Avionics  Laboratory  In  response  to  their  dosign  request,  Tho  design 
meets  all  of  the  packaging  requirements  and  exceeds  all  of  the  per- 
formance requirements  by  a significant  margin. 

A secondary  conclusion  was  that  tailless  vehicles  as  a class 
should  be  excluded  from  further  consideration  as  a possible  configur- 
ation to  meet  the  desired  design  requirements.  This  class  of  vehicles 
has  major  stability  problems  and  at  least  in  the  above  analysis,  no 
possible  solutions  to  these  problems  were  available.  It  appeared  that 
If  the  launch  dynamics  of  this  class  of  vehicles  wore  considered,  the 
stabil’ty  problem  would  oe  insolvable. 

Recommendat i ons 

It  is  recommended  that  further  analysis  of  the  proposed  vehicle 
be  conducted  Into  the  areas  excluded  from  the  above  analysis.  In  par- 
ticular, the  launch  dynamics  need  to  be  carefully  studied  because  it  is 
possible  that  the  perturbations  from  the  launch  would  cause  the  vehicle 
to  stabilize  Into  some  undesirable  flight  condition.  To  study  this  aroa, 
an  extensive  set  of  wind  tunnel  tests  needs  to  bo  conducted  to  determine 
the  vehicle  dynamic  stability  derivatives  so  that  a six  degree  of  free- 
dom non-linear  analysis  can  be  conducted.  Finally  flight  tests  of  the 
vehicle  on  a fully  Instrumented  test  range  under  simulated  operational 
conditions  need  to  be  conducted. 
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In  the  appendices  that  follow  thn  longitudinal  data  Is  normalized 
with  9*  1,0,  nd  the  lateral  data  Is  normalized  with  ♦ » 1,0  , The  dashed 
vectors  Indicate  a magnitude  too  small  to  be  visible. 
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Appendix  H 

Intern  I Reoronsn  Data  For  The  Wlnn-Oody  Conf I miration 
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Appendix  0 

Long! tudi nat  Per i vnt i vn  Sons I ti v i tv  Analysis 
Derivative:  Short  Period:  Phugold: 
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Appendix  E 

Lotera I Response  Data  for  the  Winq-Body-Tai I Conf Iquratlon 
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Appendix  F 

Latora I Dor  1 vat i ve  Sensitivity  Analysis 
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Appendix  6 

Lai  ora  I Response  for  Mortified  Wlnn-Rody-Ta  1 1 Configuration 
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* I 3738: -.0353 


Spiral  Mode 


B .0106:1 : -2. 24 18 
B:p:r  - I :-. 1061:. 2783 


Rolling  Mode 

= -.0478:1:-. 0999 
0:p:t  = I 4236:-. 0422 


Spiral  Mode 


= .0109: 1 :-2. 1634 
Bsp:r  * I:-. 1080:. 2718 


Rolling  Mode 


8:$:i|»  = -.0452:1:-.  1051 
8:p:r  = I : .44 10:-. 0462 


Spiral  Mode 


= .01 1 1: 1:-2.09I2 
= I 1098:. 2656 


gaeac , sacagra 
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Appendix  I 

Latera I Response  Data  for  the  Cruci form  Tal I Conf Ipuratlon 


Dutch  Roll  Mode 

Rol  I Ing  Mode 

■ .0018: 1 : — * 25 1 2 
6:p:f  » I t-6.3512: 1.6080 

Spiral  Mode 

* .0066: 1 :-2. 0008 
8:0:^  * I :-. 1949:. 4477 


Rolling  Mode 


0:$:C»  • .0020: I :-. 2526 
8:0:P  a I :-5. 9322: 1.5105 


Spiral  Mode 


0: ^ : ip  » .0066: 1 :- 1 .9629 
B’Ptt  * I 1988: .4465 


.30 


Rol I Ing  Mode 

9 .0021 : 1 :-.254 1 
B:p:r  « I :-5. 5077: I .4106 

Spiral  Mode 


Bt$!'}'  9 .0067; ! :-l  ,91 53 
8:{>:f  * I 2014: .4397 
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